Azimuthal asymmetries of the large-scale circulation 
in turbulent Rayleigh-Benard convection 
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Previously we published a dynamical model (E. Brown and G. Ahlers, Phys. Fluids, 20, 075101 
(2008)) for the large-scale-circulation (LSC) dynamics of Rayleigh-Benard convection in cylindrical 
containers. The model consists of a pair of stochastic ordinary differential equations, motivated 
by the Navier-Stokes equations, one each for the strength 8 and the orientation of the LSC. 
Here we extend it to cases where the rotational invariance of the system is broken by one of several 
physically relevant perturbations. As an example of this symmetry breaking we present experimental 
measurements of the LSC dynamics for a container tilted relative to gravity. In that case the model 
predicts that the buoyancy of the thermal boundary layers encourages fluid to travel along the 
steepest slope, that it locks the LSC in this direction, and that it strengthens the flow, as seen in 
experiments. The increase in LSC strength is shown to be responsible for the observed suppression 
of cessations and azimuthal fluctuations. We predict and observe that for large enough tilt angles, 
the restoring force that aligns the flow with the slope is strong enough to cause oscillations of 
the LSC around this orientation. This planar oscillation mode is different from coherent torsional 
oscillations that have been observed previously. The model was applied also to containers with 
elliptical cross-sections and predicts that the pressure due to the side walls forces the flow into a 
preferred orientation in the direction of the longest diameter. When the ellipticity is large enough, 
then oscillations around this orientation are predicted. The model shows that various azimuthal 
asymmetries will lock the LSC orientation. However, only those that act on the 5-equation (such 
as tilting relative to gravity) will enhance the LSC strength and suppress cessations and other 
azimuthal dynamics. Those that affect only the 9o equation, such as an interaction with Earth's 
Coriolis force, will align the flow but will not influence its strength and the frequency of cessations. 



I. INTRODUCTION 

The problem of Rayleigh-Benard convection (RBC) 
consists of a fluid sample heated from below and cooled 
from above (for reviews, see Refs. [H, H, Q)- In tur- 
bulent convection, a hot (cold) thermal boundary layer 
at the bottom (top) becomes unstable due to buoyancy 
and emits hot (cold) volumes of fluid known as "plumes" 
which detach from the boundary layer and transport heat 
vertically. In cylindrical containers with aspect ratio 
T = D/L 1 [L \s the height and D is the diame- 
ter of the sample) these plumes contribute to the driv- 
ing and are carried by a large-scale circulation (LSC), 
also know as the mean wind, which forms a loop with 
up-flow and down-flow on opposite sides of the sample 

0, i, S 0, S i 0, [O, d, 111 Hi, [H, M [13, p. The 

LSC breaks the rotational invariance of a cylindrical sam- 
ple and its near-vertical circulation plane must somehow 
choose an azimuthal orientation. This orientation has 
been found to undergo spontaneous diffusive meandering 
[l^ [13, [2l|, HI] . The LSC also undergoes re-orientations 
both by azimuthal rotations [3, [2l|, and by cessations in 
which it essentially slows to a stop and then restarts in a 
near-random new orientation [21, 23]. The dynamics of 
the LSC also are known to include torsional oscillations 
in which the orientation of the upper half oscillates out 



of phase with the lower half [H, [H [H, HI] . 

For cylindrical containers much of the LSC dynamics 
are described very well by a model presented in Ref. [l^] 
and studied in detail in Ref. |28i] (an exception is the tor- 
sional oscillation; although it is known that this mode has 
the properties of a damped oscillator driven by the tur- 
bulent background fluctuations 26], its origin remains ob- 
scure at this time). This model was physically motivated 
by the Navier-Stokes (NS) equations and was reduced to 
a pair of stochastic ordinary differential equations. The 
LSC strength was represented by a temperature ampli- 
tude 5. The change in 5 is determined by the balance 
between buoyancy and drag forces, and the system is 
driven by stochastic fluctuations fs{t) which phenomeno- 
logically represent the action of the small-scale turbulent 
background fluctuations on the LSC: 
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Here do is the fixed-point amplitude of the deterministic 
equation and is nearly equal to the mean value of 6 ex- 
pected to be observed in experiment. The azimuthal ac- 
celeration of the LSC orientation Oq is damped by the ro- 
tational inertia of the LSC and driven by stochastic fluc- 
tuations, again representing the turbulent fluctuations: 
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The stochastic terms were assumed to be Gaussian dis- 
tributed white noise, leading to experimentaUy measured 
diffusivities Ds and Dg. The timescales were predicted by 

the model to be ts = L'^ /{ISi^rIJ'^) and = L'^/{2vRe) 
where v is the kinematic viscosity of the fluid and Re is 
the Reynolds number. Measurements of the LSC model 
parameters from the same experimental apparatus as 
the present measurements were reported in Ref. [1^ (for 
the diffusivities and timescales) and Ref. [1^ (for the 
Reynolds number). 

The goal of the present work was to study the LSC 
dynamics in systems with slight azimuthal asymmetries 
to determine if and how results found in cylindrical con- 
tainers can be applied to such asymmetric systems. This 
is an important extension since most natural convection 
systems do not have the symmetry of ideal laboratory ex- 
periments. One such asymmetry is Earth's Coriolis force, 
which was found to cause a net rotation of the LSC ori- 
entation on average once every 3 days, and to align the 
LSC in a preferred orientation close to West [22|. One 
natural convection system with azimuthal asymmetry is 
Earth's outer core which drives Earth's magnetic field 
with orientation changes [s^ which some believe could be 
caused by cessations. Another is the solar convection 
zone which has both torsional oscillations and periodic 
reversals [3l|. Both of these examples are spherical ge- 
ometries which become asymmetric due to the action of 
external fields. Reversals are also known to occur in the 
wind direction in the atmosphere [33 |. 

In Sect, mil we extend the model of the LSC dynamics 
to the case where the sample is tilted slightly relative to 
gravity. It is known from experiments that tilting the 
sample breaks the rotational symmetry significantly and 
aligns the flow along the direction of tilt with the flow 
going up along the hot bottom plate and down along 
the cold top plate u, 33 1. It has been shown that the 
Reynolds number [3j, [SJ and average amplitude of the 
LSC [33*1 increase with tilt due to the component of buoy- 
ancy of fluid in the thermal boundary layers parallel to 
the plates, and that the average frequencies of occur- 
rence of both cessations and rotations are reduced [21 1 . 



The orientation-locking effect is due to a component of 
the same buoyancy that aligns the flow with the steepest 
slope of the plates The models of Refs. [2ll,[3l,[33 
used a steady-state force balance between buoyancy and 
drag in the boundary layers. Here these forces are gen- 
eralized to the dynamic case so they can be incorporated 
in the dynamical model of Ref. [23|. In Sect. IIVI the 
model predictions for a tilted sample are compared with 
experimental results. The good agreement provides fur- 
ther support for the model of Ref. [23]. We also report 
a prediction of a new in-phasc azimuthal oscillation due 
to the azimuthal asymmetry, and the experimental ob- 
servations of this mode, in a tilted container. Section IVl 
contains an extension of the model to non-circular con- 
tainers and comparisons with the deterministic model of 
Ref. [2^ . Section IVIBI contains a discussion of the ex- 
perimental results of Ref. [2^ on the effects of Earth's 



Coriolis force and of asymmetric heating, and how those 
results relate to the stochastic model. 

Finally we note that several other models of the LSC 
have been presented in the literature. We mentioned four 
of them in a previous paper [2^, ,3^ . [36l . [37| and discussed 
how they differ from our model; we shall not repeat this 
discussion at this point. However, very recently yet an- 
other model was developed ^] for a system with 
a square cross section of aspect ratio four with periodic 
lateral boundary conditions. It considers the coupling 
between the temperature and velocity field, but has no 
azimuthal degree of freedom. It is a deterministic model 
with relaxational dynamics, and thus can not reproduce 
any of the statistically stationary dynamics (which in- 
clude fluctuation amplitudes and cessations) of our model 
and of the physical system. 



II. THE EXPERIMENT 

Measurements were taken in the cylindrical sample de- 
scribed in detail as the medium sample of Ref. 40] . The 
top and bottom plates were made of copper and the 
side wall was made of plexiglas. The sample had height 
L — 24.76 cm and diameter D = 24.81 cm, yielding an 
aspect ratio F = 1.00. The container was surrounded by 
a thermostatic shield kept at the fluid temperature at the 
mid-height of the sample to minimize heat flux through 
the side of the sample. The working fluid was water at 
a mean temperature of 40.00°C. The Prandtl number cr, 
defined by 
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(3) 



(k is the thermal diffusivity of the fluid) was equal to 
4.38 {ly = 6.7 x 10^^ m^/s) for our experiments. The 
temperature difference AT between the top and bottom 
plates was controlled with a resistive heater imbedded in 
the bottom plate and cooled with water pumped through 
a channel in the top plate. The Rayleigh number i? is a 
control parameter defined by 
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where a is the isobaric thermal expansion coefficient and 
g is the acceleration of gravity. We ran experiments with 
different AT to obtain a Rayleigh number range of 2 x 
10® ^ R ^ xlO^". The other control parameter was the 
tilt angle /3 of the sample relative to gravity, which was 
changed by turning a leveling screw to vary the extent of 
the azimuthal asymmetry of the system. 

The LSC was observed using thermistors imbedded in 
the side wall from the outside so they were near to but 
not penetrating into the fluid region [2y, [s^ . The interior 
of the container had no measuring devices or other ob- 
jects that could break the azimuthal symmetry of the sys- 
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tern. Depending on R, the dominating azimuthal asym- 
metry for the level-sample experiments was either due to 
Earth's Coriolis force or to a small asymmetry in the top- 
plate cooling system ^2^. Eight thermsitors were placed, 
equally spaced azimuthally, at each of the three heights 
—L/A, 0, and L/4 (the origin of the vertical axis is taken 
to be at the horizontal midplane of the sample). The 
sampling period was about 2.5 seconds. Since the LSC 
carried warm fluid from the bottom plate up one side and 
cold fluid down the other side, the temperture profile in- 
dicated the orientation and strength of the LSC. The 
thermistor readings were fit, separately at each height 
and time step, by the function 



T = To + Scos [e-Bo] 



(5) 



Of the three fit parameters - the temperature amplitude 
i5, the mean temperature Tq, and the azimuthal orienta- 
tion 9o - 6 and Oq are the same two variables that were 
used in the stochastic model to describe the LSC. As de- 
fined here, the orientation Oq is on the side of the sample 
where the LSC is warm and up-flowing. Results from 
this method of determining the orientation and strength 
of the LSC were reported in several previous publications 
[m m, ia, .2.6, 28, .29, ,33, ,41, ,42, 43,]. Measurements of 
the LSC orientation and amplitude reported in this paper 
were taken at the mid-height unless otherwise noted. The 
measurements at all 3 heights are necessary to measure 
and distinguish between the torsional oscillation mode 
p"6l] and the in-phase oscillations that we found due to 
tilting the sample. 



III. THE MODEL FOR TILTED CONTAINERS 

A. The model equations 

We now consider how the LSC dynamics change when 
the sample is tilted slightly relative to gravity with the 
high side at the orientation 9p. The additional forces that 
are not present in the level sample are formulated as vol- 
ume averages so that they can be added as perturbative 
terms to the model represented by Eqs.[T]and[2l By tilt- 
ing the sample, the buoyancy of the thermal boundary 
layers is no longer aligned with the cylinder axis, but a 
component sin /3 of that force is along the direction of tilt 
Op. Fluid in the bottom boundary flows up-slope, while 
fluid in the top boundary flows down-slope, so their con- 
tributions add constructively to enhance and align the 
LSC in the direction of Bp . The additional acceleration in 
the direction parallel to the plates is given by the Navier- 
Stokes equation to be u' = ga{T — To)sin/3. Approxi- 
mately half of the temperature drop in the sample occurs 
in each boundary layer, so the average temperature of 
each of the thermal boundary layers relative to the mean 
is AT/4. The total volume fraction of the boundary lay- 
ers is 21/ L, where I = L/{2Af) is an adequate approxi- 
mation for the thermal boundary-layer width j4^] {J\f is 



the Nusselt number). Multiplying these terms gives the 
volume-averaged acceleration 



., gaAT . 
[u )v ^ — TTT' sm /? . 
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A fraction cos(0o — ^z?) of this acceleration is in the di- 
rection of the LSC and adds to the LSC strength when it 
is ahgned with 9/3, and a fraction sin(0o — ^/j) pushes in 
the azimuthal direction towards 9/}. The model describes 
the LSC strength in terms of the measured parameter S, 
so we use the assumption that 5 (x U [131 where U is the 
maximum velocity in the coordinate of the LSC (which 
occurs at the edge of the bulk region just outside of the 
boundary layers) to get 
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The azimuthal rotation rate 9o is related to the azimuthal 
velocity ug by 9o = ug/r where r is the radius in cylin- 
drical coordinates. The volume average for becomes 
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This calculation differs somewhat from those of 
Refs. [sl, [131 because in those cases the force balance 
between buoyancy and drag was made only in the viscous 
boundary layer. Following the spirit of Ref. ^27:], here vol- 
ume averages were taken, which allowed a consideration 
of buoyancy in the thermal boundary layer separate from 
the drag in the viscous boundary layer. These differences 
do not affect the qualitative dynamics of the model, but 
they do change numerical coefficients by factors of order 
1 and influence slightly the dependence of the predictions 
on the Rayleigh number. The amplitude equation Eq. [1] 
with the addition of the tilt term of Eq. [7] becomes 
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+ fs{t) + 5p sin/3cos(6io - 9p) 



(9) 



Using the definitions of R (Eq. d]), a (Eq. [3]), and the 
Reynolds number R^ = {U)L/v we find 
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The azimuthal equation Eq. [2] with the addional term of 
Eq. [5] becomes 



70 = 



S9o_ 



+ fQ{t)-ujlsmPsm{9o-9p) (11) 



where 
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(12) 



Equations [5] and [TT] compose the model of the LSC dy- 
namics in tilted containers. Most of the data presented in 
this paper are taken at i? = 2.8 x 10®. For this R, it was 
found that So = 0.020 K, ts = 133 s, Ds = 8.7 x 10"® 
K^/s, Tg = 14 s, and Dg = 1.5 x lO"*" rad^/s^ using 
the methods of Ref . [28^ . Making additional use of mea- 
surements of A/" \A5^ and ^29] from the same appa- 
ratus, the predictions for the tilt parameters at this R 
are Sfi = 1.2 x lO""* K/s and wl = 1.3 x lO""* s^^. The 



measured values of Sp and uj^ are reported in Sect. US 
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B. The potential 

The stochastic amplitude equation can be understood 
as describing the diffusion of (5 in a potential V defined 
by V = ^ J SdetdS where Sdet is the deterministic part of 
Eq. m This gives 



FIG. 1: The potential V = - J 6ddS from Eq. [13] for three 
tilt angles at i? = 1.1 x 10^° and 9o = Op. The potential is 
non-dimensionalized by ts/5q so the (3 — Q curve is the same 
for all R, but the value of Sts/Sq and hence the magnitude of 
the shift with /3 varies with R. The displacements relative to 
zero of the horizontal dotted lines indicate the heights of the 
potential barriers AV for cessations for each tilt angle. 



V = -— + - j= -SpS sin f3cos{eo- Of}) . (13) 

This potential has the form of a well with a stable fixed 
point where S{S) — 0. For f3 = this point occurs at 
6 = Sq . The potential goes to infinity for large S, but has 
the finite value V = at 5 = 0. When fitting the temper- 
ature profile Eq. O to the data, 6 is forced to be positive 
definite so that S and 6*0 can be determined uniquely. 
Thus, if d approaches zero, then this implies that the 
LSC is about to reverse direction. These reversals are 
handled with a refiective boundary condition ai S — 0, 
and upon refiection ivr is added to 9q. This is equiv- 
alent to a double-well potential in S that is symmetric 
around 5 = 0. The potential is shown in Fig. [1] for three 
tilt angles. The minimum of the potential depends quite 
strongly on (3. It is immediately clear that cessations will 
be strongly suppressed by a tilt with the deeper mini- 
mum because they depend on diffusion, driven by the 
background fluctuations, from the neighborhood of the 
potential minimum to V ^ 0. 

C. The linearized model 

For ease of calculations, the model can be linearized 
by expanding around the stable fixed point at the mini- 
mum of the potential well with the substitution e = S — Sq 
as was done in Ref. [2^. Since S fluctuates around the 
potential minimum, this is a good approximation for cal- 
culations of most average values, but it will not work for 
studying cessations which occur when 6 becomes small. 
It is known from experiments that the orientation 0q be- 
comes restricted to a narrow range of 9 near 9f3 even 



for very small tilt angles [2l[. For the smallest mea- 
sured tilt angle f3 = 0.017 rad, the orientation was con- 
fined to a range of width 0.38 rad around 9^ so that 
(cos(0o — 9p)) w 0.93. Thus, setting the cosine term 
equal to 1 underestimates the tilt term by only 7%, and 
this error diminishes rapidly for larger tilt angles where 
the orientation is confined closer to 9p. Keeping only the 
lowest-order term in e, the linearized version of Eq. [5] is 
then 



e = -— + fs{t)+6psinf]. (14) 
Zts 

This equation is valid for tilt angles large enough to lock 
the orientation but small enough that the fixed point am- 
plitude S remains in the parabolic region of the /3 = po- 
tential. The fixed-point amplitude in this approximation 
is given by 

e = 2T5Spsinp. (15) 

The linearized potential is defined by Ve{e) = V{So) — 
J edetde where £det is the deterministic part of Eq. 1141 
This gives 

V,{e) = ^- {So + e)Sp sin 13- ^ . (16) 
4:Ts IOts 

The azimuthal equation Eq. [11] can be linearized for 
small 9^ = 9o — 00. Setting S equal to the constant fixed 
point value 5 one gets 
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FIG. 2: The average amplitude (S), which is used as a measure 
of the fixed-point amplitude 5, as a function of tilt angle l3 at 
R = 2.8 X 10**. The solid line is a linear fit to the data. 
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FIG. 3: The amplitude variance cr| oi p{S) normalized by the 
zero-tilt mean-square amplitude Sq as a function of tilt angle 
l3a.tR — 2.8 X 10*. The horizontal dotted lines correspond 
to the mean values at small and large /3. 
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The fluctuations in 6 are relatively small (see Fig. so 
ignoring them will not affect the calculations of most av- 
erage values, but again the linearized equation will not 
be vahd for cessations when 5 becomes small. 



IV. RESULTS FOR TILTED CONTAINERS 
A. The amplitude 6 

The time-averaged measured amplitude {5) is shown 
in Fig. [5] for various tilt angles at i? = 2.8 x 10*. Since 5 
fluctuates in a near-parabolic well, {5) is a good approx- 
imation to the fixed-point amplitude 5. Consistent with 
Eq. [13 {S) varies linearly with (3 at these small values 
of /3 where sin /3 ~ /?. A fit of {5) = (5o + ^ sin (3 gave 
A = 0.094 K. This should be compared with the predic- 
tion oi A — 2Ts5p ~ 0.031 K from Eq.[T5l which is about 
a factor of 3 too small. An experiment with the same 
apparatus but at i? = 1.1 x 10^° found A = 0.68 K for 
small (3 dll, to be compared with the prediction Ak 0.33 
K which is about a factor of 2 too small. Considering 
the various approximations made in the derivation of the 
model, we regard the agreement within factors of 2 or 3 
as satisfactory. The predictive power of the model can 
be improved by using these fits to obtain empirical values 
5p = 0.00035 K/s at i? = 2.8 x 10** and 5p = 0.0072 K/s 
at i? = 1.1 X 10^*^. These values can be used as input for 
other calculations that depend on these parameters. 

Experimental results for the variance tT| of 5, normal- 
ized by the zero-tilt mean-square amplitude are shown 
in Fig. [3] as a function of /3 for i? = 2.8 x 10*. Since 
(j1/5q ^ 1, the assumption that the variable damping in 
Eq. HHhas little effect on average values is justified. 



Diffusion of 5 in a parabolic potential such as that given 
by Eq. (TH] with (3 — G was discussed in Ref. 28] . The 
probability distribution p{5) is a Gaussian function with 
a variance 



(jj = DsTs . 



(18) 



In the linearized approximation given by Eq. I14[ the po- 
tential is parabolic for small e. The tilt only adds a lin- 
ear term, so that the sum is a shifted parabola with the 
same curvature. Thus, the variance should continue to be 
given by Eq. [15] even for /3 > 0. The parameters rg and 
Ds on the r.h.s of that equation are system parameters 
that should not depend on the tilt angle. Thus we ex- 
pect the variance of S to be independent of (3. Consistent 
with that prediction, the experimental values are seen to 
be nearly independent of /3 for /3 < 0.19. However, they 
are seen to step down to a smaller constant for larger 
(3. The sharp drop in dX (3 ~ 0.19 rad is inconsistent 
with a continuous single-well potential. An analysis of 
thermistor time-series indicates a sudden change in the 
behavior of the two thermistors farthest from the plane 
of the LSC, at approximate orientations 0oii'/2, for tilt 
angles above /3 « 0.19 rad. For = 0, the fluctuations at 
9 « 9o±it/2 are equal to those at 9q and only about 15% 
less than the maximum which occurs at intermediate an- 
gles (see Ref. [i^ for data at /? = 0). A similar result 
is found for larger j3 < 0.19 rad. For (3 > 0.19 rad, the 
fluctuations at w ± 7r/2 become much smaller than 
the other thermistors by about 2 orders of magnitude. 
This suggests some sort of change in the LSC structure 
where the LSC becomes more steady and where perhaps 
plumes are better aligned with the LSC. 



B. Cessation suppression 

For small tilt angles, a suppression of cessations in a 
tilted sample was found in Ref. [2l|, but the rarity of 
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events allowed the collection of only a few data points. 
The rate of cessations is shown in Fig. 2] for a few small 
tilt angles at i? = 1.1 x 10"'^°, where cessations are more 
frequent than at i? = 2.8 x 10*. The error bars rep- 
resent the expected standard deviation of the mean for 
Poissonian statistics. In 12 days of measurements with 
tilt angles in the range 0.017 rad < f3 < 0.21 rad, no 
cessations were found. 

Using the model of diffusion in the potential well given 
by Eq. [T31 a cessation occurs when 6 diffuses from the 
bottom of the well up to the top of the barrier at (5 = 0. 
The potential barrier to cessations is then AV = V{0) — 
Ve{e), and our system is equivalent to the Arrhenius- 
Kramers problem. For large potential barriers, the rate 
of cessations is given by ojc (x exp{—2AV/Ds) where the 
proportionality constant depends on the curvatures of the 
potential near the peak and minimum [i^l • From Eq. [TBI 
the potential barrier in the case of small P is 



AV{(3) = Ay (0) + SoSp sin /? + tsS^ sin^ /? 



(19) 



The curvature of the potential minimum does not vary 
with /3, but the curvature near the peak does 1461 . Using 
the large-barrier-limit prediction of Kramers [47| and the 
empirical potential of Ref. 28] (in which a new constant 
_B = 12 is introduced), the frequency of cessations is given 
by 
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26/360 sin (3 
Ds 



exp 



2AV{f3) \ 
Ds J 



(20) 

The sin^ P term of AV is smaller than the sin /? term 
for sin/3 < 0.75, and the sin/3 term in the coefficient is 
smaller than the constant term for sin/3 < 0.20, so for 
(3 ^ 0.20 rad Eq. [20] simplifies to an exponential decay 
with p. 

For /3 = the predicted rate of cessations agreed with 
the experimental result [28] within a factor of order one. 
However, in order to better test the /3-dependence of 
the model, the proportionality Wc oc exp(— /?//5o) with 
/3o = Ds/ {260613) was fit to the data, adjusting the pro- 
portionality constant while fixing /3o = 0.017. This value 
of Po is obtained using the prediction of 6/3 from Eq. [TO] 
and the experimental values for Ds and <5o that were ob- 
tained in Ref. [281 for i? = 1.1 x lO^*'. The result is 
the solid line in the figure. The approximate agreement 
with the data indicates that the exponential decay rate is 
consistent (within a factor of two or so) with the exper- 
imental results. Adjusting the proportionality constant 
as well as Pq gave Po = 0.010 and the dashed line in the 
figure. The value is not very far from the value 0.017 
expected on the basis of the independent determinations 
of the parameters. Although the poor statistics of the 
frequency of cessations does not permit a very stringent 
test of the model, it is noteworthy that both experiment 
and model indicate that even a very small tilt will almost 
completely suppress the cessations. 
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FIG. 4: The rate of cessations ujc as a function of tilt angle 
at i? = 1.1 x W. Solid line: A fit of loc oc exp(-25o5/3/3/-D4) 
to the data with Dg / {2S0S13) fixed at the value 0.017 based 
on measurements at /? = 0. Dashed line: a similar fit, with 
Dg /{25oSp) — 0.010 least-squared adjusted. 



C. Azimuthal fluctuations 

The stationary probability distribution p{0^,) can be 
found from the steady-state Fokker-Planck equation cor- 
responding to Eq. [17] which balances the advection and 
diffusion of probability. In Eq. 1171 in addition to the ex- 
plicitly /3-dependent tilt term, the damping term depends 
on P because 6 is /^-dependent. The tilt term will be as- 
sumed to be small compared to the damping term so it 
can be ignored. This is valid in the limit ^ tOpPag 
which is satisfied by the data. The resulting Fokker- 
Planck equation gives 



To6( 
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(21) 



for the statistically stationary state. 
p{9*) oc exp(— ^J/(2(t|)) with variance 



The solution is 



(22) 



The peak of p(6'* ) was found to be Gaussian and for /3 = 
the predicted variance tr? was found to agree remarkably 

well with experiment j28|. The value (t| is plotted as 
solid circles in Fig. [S]for varous tilt angles at i? = 2.8 x 
10*. It is seen to decrease with /3. The /3-dependence 
of the predicted value is expected to be contained in (5, 
so the value cr'^^6/60 is expected to be independent of /3. 
However, the open squares in Fig. [5] show that this value 
also decreases with P, albeit more slowly. Multiplying 
cr| by {5/ 60)^ J a best fit to a constant gives n — 2.3, 
but we have no explanation for the large error in the 
predicted exponent. While this represents a significant 
disagreement of the model with the data, this width does 
not affect the overall phenomenology of the model. 
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FIG. 5: The variance of rotation rate (t| x (J/^o)" as a func- 
tion of tilt angles /3 at R — 2.8 x 10*. The variance is mul- 
tiplied by (S/So)" with several different values of n. Solid 
circles: n = 0, i.e. the bare variance (t|. Open squares: 
n = 1, which is predicted by the model in the linear approxi- 
mation to yield a constant independent of /3. Open triangles: 
n = 2.29, where n is obtained from a best fit of a horizontal 
line to the data. 



The power-law relationship between cr? and S is more 
general and not restricted to tilted samples. The Gaus- 
sian probability distribution and variance from Eq. [221 
apply instantaneously (with 5 replaced by 5) if the fluc- 
tuations in ^0 are strongly damped so that a station- 
ary p{9o) is reached in a time short compared to that 
over which S changes. This is true if Tg <^ Tg, which is 
satisfied for the experimental data. The near-Gaussian 
p{()o) implies |^o| — ■\/2/7r cTg; so it is expected that 
l^ol oc S"'^^'^. The measured average relationship was 
found to be |^o| oc (5-i i6±o.06 fj-Qjj^ series at /3 = 
for a wide range of R [2lj . This gives the exponent for the 
variance 2.32 ± 0.12 which is consistent with the scaling 
relationship in the tilted sample. The twist displacement- 
angle was found to have a similar power-law relationship 
|0'((5)| cx which suggests a displacement 9' may be 

damped by the same mechanism as Oq. 

The same scaling for with d as the instantaneous 
relationship suggests that the reduction in ag with /3 is 
indeed due to the increase in 6 as predicted. In other 
words the suppression of azimuthal fluctuations is due to 
the increase in inertial damping with 5 and is unrelated 
to the orientation-locking term of Eq. [TlJ The reduction 
in ffg due to the increased amplitude forcing must be the 
cause for the suppression of reorientations for small tilt 
angles that was found in Rcfs. [2l], [s^. The azimuthal 
forcing was measured and found not to be strong enough 
to suppress relatively fast reorientations for small tilt an- 
gles despite being able to lock the LSC orientation [2l| . 
For a large enough tilt, the orientation- locking of the az- 
imuthal equation may play a larger role in suppressing 
reorientations. 



D. Oscillations in a tilted sample 

The linearized azimuthal equation Eq. [17] is that of a 
stochastically forced damped harmonic oscillator. The 
power-spectral density P of its solutions is given by 



P 



[Su;/{6oTg)]^ + {l 



i sin /3 - 



(23) 



It has a maximum at the resonant frequency w,- = 
ujp sin P — [S/^doTg)]"^ /2 provided that ujr is real. This 

will be the case when sin/3 > [S / {SoTgLUfs)^ /2. This un- 
derdamped regime was not attainable in the parameter 
range of the experiment because the smallest measured 
value of the damping parameter was l/r^ « 0.070 s^^ at 
R — 2.8 X 10^. This is still large compared to the restoring 
frequency ajfj ~ 0.011 s~^ even for the maximum possible 
tilt P = tt/2. 

Looking back at the full azimuthal equation, Eq. 1111 
one sees that the variable damping term (containing S{t) 
rather than S) can add signiflcantly to the excitation of 
oscillations when the damping occasionally fluctuates to 
smaller values so that the restoring term temporarily 
overcomes the damping. This feature can be incorpo- 
rated in Eq. [iTj by adding the variable damping in the 
form of a stochastic term: 



(24) 



Near the stable fixed point one sees from the linearized 
amplitude equation Eq. [TJj that the multiplicative noise 
is ri{t) — e{t)/So. This means r]{t) is Gaussian distributed 
with a variance — (Jg/Sg. It is not white noise but has 
the same correlation time = 2ts as e(i) p8l |. In other 
words, this is colored noise of the Ornstein-Uhlenbeck 
form with correlation function 



(r?(<)77(t + r)) =a2cxp 



(25) 



The power spectrum corresponding to Eq. [Mj can be 
described by a complicated expression that will not be 
reproduced here [48|. A more physically clear result 
is obtained in the case where ri{t) represents Gaussian 
white noise. In this case the term in the denominator 
of Eq. [221 due to the damping is reduced by a factor of 
1 — cTgTs / {6QTg) so that the resonance condition is loos- 
ened and peaks become stronger with the fluctuation 
strength as- The value of crj^s / {SoTg) w 0.5 depend- 
ing on R. This is in a range where factor of two dif- 
ference in the fluctuation strength allows the system to 
go from overdamped to effectively undamped, and given 
the approximate nature of our calculation we cannot with 
any confidence predict values of the power spectrum or if 
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FIG. 6: Time series of the orientations at the 3 heights 
at R = 2.8 X 10^ (a): for /? = 0.09 rad. (b): for /? = 0.26 
rad. Dashed line: 6o{L/A). Sohd line: Oo{0). Dotted line: 

eo{-L/-i). 



there will be a resonance peak. This conclusion applies 
when we use the full calculation of the power spectrum 
for colored noise as well. We can only rely on the pre- 
diction that in-phase oscillations are likely to found only 
at large tilt angles with a resonant frequency that ap- 
proaches the natural frequency LOp^/sm/S from below as 
the resonance peak becomes sharper. 

We searched for these azimuthal oscillations experi- 
mentally at various tilt angles up to /? = 0.26 rad at 
i? = 2.8xl0^. A complicating factor is that an un- 
related torsional oscillation exists for turbulent convec- 
tion in cylindrical containers even for /3 = 0, in which 
the orientations of the upper and lower portions of the 
LSC oscillate out-of-phase around the mid-plane orien- 
tation [l^, [iHl- The oscillation predicted due to tilt is 
expected to have a ^-independent phase (except possibly 
for some random fluctuations) , that is the orientations of 
the lower, middle, and upper portions of the LSC are 
all aligned while they oscillate around the orientation 
Op. Since the torsional oscillation is anti-symmetric in 
z, the signal of 9^, = 9o — 9p at the mid- height should 
only have oscillatory components if the new in-phase os- 
cillatory mode exists. The out-of-phase oscillation is ob- 
tained from differences between the orientations at dif- 
ferent heights, such as the signal 9' = 9o{L/'i) — 9o{0). 

Time series of the orientations 9o{z) at the 3 heights 
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FIG. 7: Correlation functions of various orientation differ- 
ences at i? = 2.8 X 10^". Solid lines: autocorrelation (?*(r) 
of Of. Dashed lines: auto-correlation (?'(r) of 6'. The sig- 
nal corresponding to the lower-plane orientation relative to 
the mid-plane orientation is similar to 6' so is not shown for 
clarity. Dotted lines: cross-correlation gtbif) corresponding 
to the torsional oscillation, (a): {3 = 0.09 rad. (b): (3 = 0.26 
rad. 



are shown in Fig. [S^ for (3 — 0.09 rad and Fig. for 
13 = 0.26 rad, both at i? = 2.8 x 10^. For (3 = 0.09, the 
top- and bottom-row orientations oscillate out-of-phase 
around the middle-row orientation, so this is the torsional 
oscillation that was observed previously. For /3 = 0.26, all 
3 rows oscillate in-phase, and at a much lower frequency 
than the torsional oscillation in Fig. [S^. 

Autocorrelations (7*(t) of 9^. and ^'(t) of 9' for j3 — 0.09 
rad and (3 0.26 rad at i? = 2.8 x 10^" are shown 
in Fig. [7^ and b, respectively. Also shown is a cross- 
correlation (r) = {[9t{t) - 9m{t)][9i{t + t) -9m[t + 
T)])/ggi{fi), which should only oscillate in cases where 
the LSC undergoes torsional oscillations. For (3 = 0.09 
rad, both ge'ir) and gtb{T) oscillate as expected for tor- 
sional oscillations. The minimum of gtb{T) at r = indi- 
cates that the top and bottom rows are out-of-phase, thus 
showing that this mode corresponds to a torsional oscil- 
lation. The lack of peaks in gg^ indicates that there is no 
discernible in-phase oscillation at (3 = 0.09 rad. This plot 
is qualitatively the same as the results of Ref. [1^ for an 
untilted sample. For /? = 0.26 rad, the cross-correlation 
(7tfc(r) is small everywhere, indicating that the torsional 
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FIG. 8: Power spectral densities. Solid lines: the Fourier 
transform of gtbir) corresponding to the torsional oscillation. 
Dotted lines: The power spectrum of 0* corresponding to the 
in-phase oscillation. Upper data: (3 — 0.09 rad. Lower data: 
/3 = 0.26 rad. Dashed line: fit of Eql^Hlto the spectrum for 
the planar oscillation at /3 = 0.26 rad. 



oscillation has nearly disappeared. Since the physical ori- 
gin of the torsional oscillation is unknown, it is difficult 
to say why this mode is suppressed by tilting the sam- 
ple. The new long-period oscillatory mode is apparent in 
both g'{T) and (7*(r). 

Power spectra of 0* for /3 — 0.09 rad and /3 — 0.26 rad 
at i? = 2.8 X lO""^" are shown in Fig. [51 Also shown are the 
power spectra corresponding to the torsional oscillation 
given by the Fourier transform of gtbir). The peak of 
each power spectrum was fit by a Lorentzian function 



1 + {lu- uJoY/al 



(26) 



with an empirically chosen exponential background for 
the in-phase oscillation and power-law background for 
the torsional oscillation. 

Since the analytic form of the power spectrum and any 
expression for its peak frequency are complicated in the 
case of the in-phase oscillations, we compare wq to the 
model prediction for the natural frequency LOp^/svaf} even 
though the resonant frequency tends to be shifted to a 
somewhat lower value than the natural frequency. The 
agreement between the resonant frequency and natural 
frequency is expected to be better when the resonant 
peak is sharper, which is found to occur for larger (3 and 
smaller R. A plot of the fit values as a function of 
tilt angle is shown in Fig. [9}i for _R = 2.8 x 10^. A fit by 
the equation = i^Jp^/SaTp is shown as the solid line. 
The expected \/sin /3 dependence is reasonably consistent 
with the data. The coefficient was found to be ujp — 
1.9 X 10~^ s~^, within a factor of two of the predicted 
value ijjfj ~ 1.1 X 10^^ s^^. Since the resonance peaks are 
small at small /3, they are likely in the range where the 
resonant frequency is shifted below the natural frequency, 
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FIG. 9: The resonant frequency ljo obtained from the peak 
of the power spectrum of 6', on log- log scales, (a): Data as 
a function of tilt angle /3 at = 2.8 x 10*. Solid line: A fit 
of cjo = cj/3\/sin 13 to the data. Dotted line: a best power law 
fit allowing the exponent to vary gives an exponent of 0.69. 
Dashed line: a best power law fit to data for /3 > 0.19 rad 
allowing the exponent to vary gives an exponent of 0.30. (b): 
Data as a function of 7? at /3 = 0.26 rad. Solid line: A power- 
law fit with the exponent fixed at the predicted value 0.355. 
Dotted line: a best power law fit allowing the exponent to 
vary gives an exponent of 0.41. 



so it might be appropriate to only fit f the data for large 
(3 > 0.19 rad which gives an exponent of 0.30. On the 
other hand, a best power law fit to all of the data gives 
an exponent of 0.69, so the data are consistent with a 
wide range of power laws between those two values. 

A plot of ujQ for various i? at /3 = 0.26 rad is shown 
in Fig. [HId. Over the i?-range of the data the Nusselt 
number can be described well by A/" = 0.159i?°-^''° [HI. 
Thus from Eq. [12] one finds luo oc i?" '^^^. A power-law 
fit with this exponent gave the solid line in the figure, 
which is in excellent agreement at smaller R where the 
resonance peaks are sharpest. The fit gave a prefactor of 
9.8 X 10~^ s~^ which is again within a factor of two of 
the predicted value 5.7 x 10~^ s~^. A best power law fit 
to the data gives an exponent of 0.41, quite close to the 
model prediction. Since the measured frequency of in- 
phase oscillations is within a factor of two of the predicted 
value with scalings in P and R that are consistent with 
the model, it is clear that the observed oscillations are 
due to the restoring force created by tilting the sample 
as predicted. 

The mean-square amplitude of the oscillations is equal 
to the power given by the integral over the peak of the 
power spectrum, i.e. to J^^^^^ P{uj)duj = irPoap. This 
power is plotted in Fig.[Tn|for both the torsional and the 
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FIG. 10: The power vrPofp equal to the integral of the peak 
of the power spectrum corresponding to the azimuthal oscilla- 
tions. Open circles: planar oscillation. Solid circles: torsional 
oscillation. 



in-phase oscillations. The power of the torsional oscilla- 
tion drops sharply at /3 ~ 0.19. At larger /3 the planar 
oscillation has the larger amplitude. This transition was 
clearly visible in time series of 9q, where only the tor- 
sional oscillation was visible for /? < 0.19, while only the 
in-phase oscillation was visible for larger f3. This transi- 
tion occurs at the same /3 as the drop in cr^ (see Fig. [3]), 
further suggesting that there may be a structural change 
of the LSC at this tilt angle. 



V. PRESSURE FORCING IN CYLINDERS 
WITH NON-CIRCULAR CROSS SECTION 

Containers with strong azimuthal asymmetry are 
known to align the LSC plane parallel to the longest di- 
ameter [IS, Ea] , but other effects on the dynamics have 
not been well-studied. Here we model the LSC in a con- 
tainer with a non-circular horizontal cross section, corre- 
sponding to an azimuthally varying diameter D(9). This 
model predicts a restoring force that is a consequence of 
pressure gradients due to the non-circular cross section 
of the sample which vanish in the circular limit. This 
restoring force leads to a preferred orientation of the LSC 
aligned with the longest diameter(s) of the container. If 
the azimuthal asymmetry is large enough the model pre- 
dicts the existence of oscillations that are in phase over 
the entire height of the sample. 

The methods of Ref. [28| are used to produce estimates 
for volume-averaged terms that can be added as a per- 
turbation to the Langevin equation Eq. [21 As usual, the 
model assumes a singe-roll LSC. It is assumed that the 
azimuthal average of D{6) is equal to L and that D{9) 
remains close to L at all so that the LSC has an as- 
pect ratio approximately equal to one, independent of 
Oq. We will only consider the flow in a central band near 
the LSC orientation where 9 ~ 9q. This simplifies the 
problem, and it is expected to be the region where the 



side-wall-generated pressure gradient has the largest in- 
fluence. The velocity in the plane of the LSC in this 
region is assumed to be 



m 



(27) 



where is the radius in the plane of the LSC and ijJ^ is 
the angular turnover frequency. For the circular cross sec- 
tion the radi al p rofile is supported by the measurements 
of Refs. [O, [l2i- The factor L/D makes the velocity at 
the edge of the bulk the same at all 9, which is the case 
if the viscous boundary-layer width and profile are inde- 
pendent of 9o- This also makes lo^ uniform. The LSC is 
forced into this nearly circular loop due to the pressure 
from the side walls, which requires the pressure term from 
the NS equation to provide the centripetal acceleration 



Ur 



dV 



(28) 



where V is the pressure per unit density. Using Eq. [27] 
this can be integrated to get the pressure field 



V 



2D{9Y 



(29) 



Now the ^-derivative is taken to get the azimuthal accel- 
eration: 



ue 



IdV _ u^lL^rlD'{e) 
' r d0 D{9fr 



(30) 



where D'{9) = dD/d9 can have a significant fractional 
variation with 9 despite D{9) remaining close to L. The 
volume-averaged azimuthal acceleration 6g which can be 
added to the Langevin equation Eq. [2] can be obtained 
by canceling out r^/r to get 



D 



(31) 



When added as a perturbation to Eq. [51 the geometry 
effect can be expressed in terms of diffusion in a potential 
given by 



(32) 



There are two fixed points where D'{9q) = 0. The longest 
diameter corresponds to the lowest potential, and thus 
its orientation is the preferred one for the LSC circu- 
lation plane. Qualitatively this is similar to the ex- 
perimental observations for samples with a square or 
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rectangular cross section [2J, |40| where the alignment 
is along the longest diagonal. Since D{9) must be tt- 
periodic, there is at least one longest diameter in the 
cell and an even number of opposite preferred orienta- 
tions along these diameters. Containers with rectangular 
cross-sections have four preferred orientations. It is in- 
teresting to note that transitions between opposing pre- 
ferred orientations would correspond to the reversals that 
have often been discussed in the literature [l3,[H,[33| but 
that are not prevalent in the circular cross-section limit 
where azimuthal diffusion in the absence of potential ex- 
trema leads to equal probability for all LSC orientations. 

The forcing came out of the pressure term as a result 
of the assumed azimuthal variation of the internal flow 
field. If on the other hand, is independent of and uj^ 
varies with 9, then there is no azimuthal force. However, 
for such a profile to satisfy no-slip boundary conditions in 
a non-cylindrical container would require the boundary- 
layer profile to vary with 9. This variation would have 
to be drastic if the variation in D{9) is comparable to or 
larger than the boundary-layer width. 

An alternate derivation of the azimuthal forcing that 
does not depend on the internal velocity profile is ob- 
tained using a simple physical model of a vector force- 
balance at the wall. The pressure at the side wall 
provides an acceleration that can be broken down into 
two orthogonal components: the centripetal acceleration 
iir = 2u^/D, and the azimuthal acceleration iig. Since 
the pressure is normal to the wall surface, simple geome- 
try gives the ratio of these acceleration comonents to be 
D{9o)/D'i9o), so 



70 = 



2ue _ ^ulD'{9o) 
D D^{9o) 



(33) 



In this description the result of a forcing proportional 
to D'{9) is not dependent on an assumption about the 
internal flow field or the assumption of small deviation 
from azimuthal invariance. However it still depends on 
M0, and is equal to the previous result using Eq. to 
get ~ at the wall. In this description it is 

clear that the pressure of the side wall - whose gradient 
is no longer directed towards the sample center for a non- 
circular cross-section - provides the acceleration to align 
the flow with the longest diameter. 

The probability distribution of the orientation can be 
obtained from the steady-state Fokker-Planck equation 
in the strong-damping limit and is 



p{9o) (x exp 



(34) 



The rate of switchings between different preferred orien- 
tations for various cross sections can be found by directly 
applying Kramers' prediction. 



A. Oscillations in containers with an elliptic cross 
section 

Let us take for example a sidewall shape that is slightly 
elliptical so the diameter varies as D{9)/ L = l-|-ecos(20) 
where e <C 1 so the aspect ratio remains close to 1. Then 
D'{9)/L — — 2esin(2f?o) where 9a = Q and tt are the pre- 
ferred orientations and 2e is the difference between the 
minor and major axes. Taking the small-angle approxi- 
mation sin(26'o) ~ 2^0 and adding this forcing due to the 
elliptical sidewall to the linear equation Eq. [T7] gives the 
stochastically driven linear harmonic-oscillator equation 



Aeul9a + fg{t) 



(35) 



Here the ellipticity parameter e plays a similar role to 
the tilt angle (3 in Eq. [T71 The qualitative difference 
between the two cases is that the container geometry is 
not expected to affect the LSC strength because it does 
not enter into the 5 equation. The resonant frequency of 
the deterministic system is 



1 

2tI 



(36) 



so the system becomes underdamped above a critical el- 
lipticity Sc = 1/ (Sui'^T^). For i? = 1.1 X 10^", for instance, 
we find Ec — 0.16. But as for the case for the tilted sam- 
ple, the inclusion of the factor 6/ So in the damping term 
of Eq [35] will lead to oscillations for e < €c- A detailed 
consideration would lead to a result similar to Eq. [24l 

The linear oscillator equation applies when the orien- 
tation is confined to a small region due to the asymmetry. 
If the asymmetry is weak, so that the potential barrier 
AVg equal to the difference between the maximum and 
minimum of the potential is not large compared to the 
effective fluctuation energy DgTg/2, then the orientation 
can switch between potential minima. If the fluctuation 
energy is much larger, switchings become frequent and 
the effect of the potential is minimal. For the elliptical 
container, the well depth is AVg — 2ea;^/(l — e^)^, so to 
have bounded oscillations, the ellipticity must be greater 
than Er « T^Dg/{2ujl) = 2r|Dgec. For R = 10^° for 
instance this yields Sr « 0.02ec ~ 0.003. Using the mea- 
sured values from Ref. [l^ one finds a very weak depen- 
dence on R. This small but non-zero ellipticity required 
to lock the orientation is consistent with experiments in 
non-circular containers which found a locked orientation, 
as well with measurements for nominally circular contain- 
ers - which must have some small effective ellipticity - 
that found that the orientation meanders. When e ^ Sr, 
then the rate of reversals between the two orientations 
along the longest diameter is given by the Kramers re- 
sult 



OJr 



2V2e 



■ exp 



(37) 
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B. Containers with rectangular cross section 

Common experimental geometries besides circular 
cylinders have been samples with square or rectangular 
cross-sections. Although the model should work best for 
smooth cross-sections with small D'{6)/L, it is instruc- 
tive to apply it to the rectangular case. With a square 
cross-section, for example, the shortest and longest di- 
ameters can be compared to get an effective ellipticity 
Eeff « (%/2- 1)/{V2 + 1) « 0.29. This is larger than the 
predicted Ec required to get oscillations. Because of the 
approximate nature of the prediction, at best we can say 
that there is a good chance of finding azimuthal oscilla- 
tions in a square cell due to the pressure of the side wall, 
and that reversals should be extremely rare or absent. 

The Fokker-Planck analysis on rectangular cross- 
section containers leads to an unusual result for p{0q). 
Consider a rectangular container of side lengths a+ and 
a_ so that the longest diameter is = (a^ -I- a^)^/^ 
with cross-section aspect ratio A = a+/a^. Let 9 = 
along a diameter and D± correspond to a wall with length 
a±. Then the diameter can be expressed as 



D±(e) 



cost 



+ A±i 



smf 



(38) 



provided \9±\ < tt— 2 arctan (A'Fi^ . This result applies to 
all four corners but is mirror-imaged at alternate corners. 
Examples of the full potential Vg given by Eqs. [32] and 
1551 are shown in Fig. [11] 

The variation in diameter is at least as large as that 
of the square cross-section, so reversals will usually be 
strongly suppressed in a rectangular container. How- 
ever, switchings between nearby corners may occur for 
sufficiently large A. For A close to one, a small-angle ex- 
pansion of p{Oq) around the longest diameters is a good 
approximation. A first-order expansion of the diameter 
around a corner gives 

D±{eo) - {1 ~ A^^\e\)D^ . (39) 

The steady-state Fokker-Planck equation Eq. [34] gives 

/ 2|0o|A±i^2^2\ 



P±{()o) « exp 



(40) 



This distribution has a width of 2er£'^/(L^A^^) which 
is on the order of 1 degree. One sees that a sample with 
a square cross section (A = 1) gives an exponential peak 
near a corner with equal decay rates on the two sides ad- 
jacent to the corner. A rectangular cross-section (A > 1) 
also gives an exponential peak, but with a faster decay 
rate on the longer side. In comparison, an elliptical cross 
section gives a Gaussian peak due to the parabolic shape 
around the longest diameter. Thus a test of this model 
would be to measure p{Oq) in containers with rectangular 
cross-sections. A polygon with corners having more ob- 
tuse angles would be expected to give wider distributions 
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FIG. 11: The azimuthal potentials Vg for containers with 
rectangular cross-sections. Solid line: A = 1, i. e. a square 
cross-section. Dotted line: A = 2. 



of p{9) that might be easier to resolve, and increase the 
rate of switchings between nearby corners. 

To determine the rate of switching between adjacent 
corners, we compare the potential barrier AVs between 
the corners with the diffusive energy TgDg/2. The poten- 
tial barrier between adjacent corners is 



2i?2 A±2 



(41) 



where A^ is used for the barrier between the closer cor- 
ners and A~2 is used for the barrier between the further 
corners. The potential barrier is equal to the diffusive 
energy when A « L/(2D„^) « 6.3 for R = 10^° 
and assuming L = D,n/y/2. Thus for larger A the 
rate of switchings between adjacent corners becomes fre- 
quent, while for smaller A the rate of switchings drops 
off rapidly with the usual exponential form proportional 
to eyi\)[—2/Ws/{TgDg)\. For large A, the orientation 
may not be locked into a single corner and the rate of 
switching between nearby corners could become appre- 
ciable. However, such large A values are far from a nearly 
uniform-diameter cross-section; thus it is unclear how re- 
liable our approximations are in this limit. 



C. Comparison to prediction for an ellipsoidal 
container 

The deterministic model of the LSC presented in 
Ref. [m used the Boussinesq equations with slip bound- 
ary conditions for ellipsoidal samples. Convection was 
due to an applied uniform vertical temperature gradient 
along one axis of the ellipsoid. It is apparent that this ge- 
ometry is not very closely related to any that have been 
studied experimentally so far. The main virtue of the 
model is that it allowed an exact solution for the velocity 
field, but dissipation had to be introduced a posteriori 
and phenomenologically. For certain parameter ranges 
the model leads to an in-phase azimuthal oscillation for 
£ > £c — 0.02 of the LSC around a preferred orienta- 
tion in which its circulation plane is aligned on average 
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along the longest horizontal diameter of the container. 
The value Ec = 0.02 is below our stochastic-model pre- 
diction for underdamped oscillations. For e > 0.04, the 
deterministic model also found various chaotic modes. A 
common feature of our stochastic model and the deter- 
ministic one of Ref. [1^ is that in both cases the oscilla- 
tions are a consequence of pressure gradients due to the 
non-circular cross section of the sample. However, in our 
model they are those of a damped oscillator driven by a 
stochastic force. 

The experimental results of Ref. [2^] included an az- 
imuthal oscillation mode with a frequency that matched 
the turnover frequency of the LSC. The measurement 
was made only near the top plate. Thus it could not 
be determined whether there was a phase difference be- 
tween the top and bottom of the flow. This was reported 
to be consistent with the in-phase oscillation predicted 
by the deterministic model presented in the same paper, 
provided that e was used as a tunable parameter. It 
seems more likely that the observed oscillations were the 
torsional oscillation mode reported in Ref. [l6j , which is 
already known to have the same frequency without the 
need for a tunable parameter. Physically, there was no 
reason to assume the large ellipticity required to obtain 
oscillations synchronous with the turnover frequency be- 
cause the container had a nominally circular cross sec- 
tion. The torsional oscillation is known to occur in az- 
imuthally symmetric containers, and because of the sym- 
metry, it does not oscillate around a fixed orientation but 
rather around the mid-plane orientation which meanders 
[11]. The experiments of Refs. [H] and [H] found that 
the azimuthal oscillations occurred around a fixed ori- 
entation in the experiment, but that may well be due 
to locking of the LSC orientation due to various small 
asymmetries which occurs to some extent in all samples 



VI. OTHER PERTURBATIONS 

A. Asymmetric heating 

It was found in Ref. [1^ that asymmetric heating can 
align the LSC orientation. A small asymmetric heating 
effect was measured in Ref. 



forcing in the direction of 9c in terms of acceleration is 



22| and found to produce 



alignment of the LSC in the direction of a small horizon- 
tal temperature gradient at the top plate. Asymmetric 
heating enhances the buoyancy of the boundary layers 
near the hotter side of the sample, and the dynamic ef- 
fects on the LSC are expected to be similar to the effect 
of tilting the sample. Thus one expects the potential 
well of 5 to become deeper, leading to a suppression of 
cessations. 

Consider a temperature difference 6T over a length L 
horizontally across a plate such that the hot side is in 
the direction Oc so the LSC is forced in this direction. 
We volume-average forces using the same methods as 
Ref. |27[ to obtain the effects on both 9q and 5. The 



{u')v = 



gaST 
37r7V 



One vector component gives the azimuthal forcing 



(r> 
ga6T 



smiOc - Oo ) 



(42) 

(43) 
(44) 



which acts to align the fiow and can be added as a per- 
turbative term to Eq. [2] The other component acts to 
strengthen the LSC when it is aligned with the asymme- 
try and gives 



= TZ C0S( 



u 

gaSTLS 



cos{6c - Oq) 



(45) 
(46) 



which can be added as a perturbative term to Eq. [1] 

B. The Coriolis force 

The experiment of Ref. [2^] studied the effects of 
Earth's Coriolis force on the LSC. The model of the 
Coriolis force presented in Ref. [23| was equivalent to 
the strong-damping limit of Eq.[2]with additional forcing 
terms that both align the LSC and cause a net azimuthal 
rotation of the LSC plane. A more complete analysis 
of the integration of the two models was presented in 
Ref. 

A significant qualitative difference between the several 
asymmetries that we have discussed is that the Coriolis 
force and non-circular side walls only add terms to the 
azimuthal equation Eq. [2] while the heating and tilt asym- 
metries add terms to both the azimuthal equation Eq. [2] 
and the amplitude equation Eq. [TJ In general the Cori- 
olis force always acts perpendicular to a flow direction. 
Thus it can deflect a flow but not strengthen or weaken it. 
Since the Coriolis force only affects the azimuthal motion 
and not the LSC strength, it should not affect cessations 
which we found to be suppressed by the increase of the 
amplitude in the case of a tilted sample. This was studied 
in Ref. [11] by tilting the container slightly and observ- 
ing that the peaks or minima of measured parameters did 
not occur at /? = 0, but at a slightly shifted angle due to 
the interplay between Coriolis and tilt forces in one set of 
experiments, and the interplay between a heating asym- 
metry, Coriolis, and tilt forces in another experiment. 
The conclusion was that the measured azimuthal forcing 
as a function of smd the parameters corresponding to 
p{0o) were found to be affected by all three mechanisms. 
However, S, Re, cTg, and the frequency of reorientations 
were only affected by the heating and tilt asymmetries. 
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This is in agreement with the stochastic model, since S, 
Rf,, CTg, and the frequency of reorientations should only 
be affected by an asymetry that also enhances the LSC 
strength. The data of Ref. [l^l were taken for only very 
weak asymmetries (the apparatus was stationary on the 
surface of the Earth, so the Coriolis force measured was 
due to Earth's rotation), so those data support the con- 
clusion that a weak Coriolis force does not suppress reori- 
entations. This conclusion has not been tested in systems 
with stronger asymmetries. A large enough asymmetry 
could change the structure of the LSC, in which case the 
models presented here assuming a single-roll LSC would 
not apply. 



VII. SUMMARY AND CONCLUSIONS 

In a previous paper [28| we developed a model for 
the large-scale circulation in turbulent Rayleigh-Benard 
convection in cylindrical containers. It consisted of two 
stochastic ordinary differential equations, one for a mea- 
sure 6 of the strength (Eq. [T|) and the other for the az- 
imuthal orientation 9q (Eq. ^ of the flow. This model 
preserved the rotational invariance of the ideal cylindrical 
system. In many, if not most, real physical systems, in- 
cluding idealized laboratory experiments, this invariance 
is broken by a variety of perturbations. In the present 
paper we considered the effects of tilting the cylinder 
axis relative to gravity, deviation from a circular cylinder 
cross section, asymmetric heating of the top plate, and 
the influence of the Coriolis force due to rotation. We 
did this by adding physically motivated terms to either 
or both of the two equations as required by the nature 
of the perturbation. Predictions based on this perturbed 
model were compared with experiments for tilted sam- 
ples, and generally there was agreement within factors of 
order unity between the measured and predicted quanti- 
ties. Predictions for the influence of the Coriolis force due 
to Earth's rotation had been compared with experiment 
in a previous paper [22|. 

A particularly interesting result of this work is that 
asymmetries can be categorized based on how they affect 
the flow dynamics. Those that influence the buoyancy of 
the thermal boundary layers, such as a tilt or asymmetric 
heating, increase the strength of the LSC. They lead to 
an additive term of the (5-equation. We show that they 
suppress cessations and azimuthal fluctuations because 
they deepen the potential well corresponding to the S- 
equation. Other asymmetries can result in the addition 
of a perturbative term only to the equation. Examples 
include the influence of a Coriolis force due to rotation of 
the sample, and the influence of deviations from a circular 
cross section of the sample. These asymmetries align the 
LSC, but do not affect its strength. As a consequence, for 
instance, they do not change the frequency of cessations 
of the LSC. 

Another interesting feature of the model is that it pre- 
dicts the generation of an azimuthal oscillation of the 



LSC circulation plane when the tilt is sufficiently strong. 
These oscillations are due to an azimuthal restoring force 
that arises when the sample LSC orientation deviates 
from the tilt direction. This case was investigated in de- 
tail by experiment, and the predicted oscillations, with 
frequencies close to the predicted values, were indeed ob- 
served. It is important to note that this oscillation is 
expected to be in phase over the entire sample length, in 
contrast to a quite different torsional mode of unknown 
orig in that is known to exist even for the untilted sample 
[2a |. Whereas the torsional mode has a period that is 
synchronous with the LSC turnover time, this new mode 
generally has a longer period that depends on the tilt 
angle and is not related to the turnover time. 

An oscillatory mode, coherent over the entire sample 
height, is predicted as well for a sufficiently large eccen- 
tricity of an elliptical cylinder cross section. In this case 
the oscillations are predicted to be due to a pressure gra- 
dient induced by the deviation from circular symmetry. 
Such an asymmetry is predicted to lead to two preferred 
orientations, with the LSC circulation plane aligned with 
the longer axis of the ellipse. Oscillations can occur about 
one or the other of these alignments. When the eccen- 
tricity is small, the LSC orientation is expected to diffuse 
in the neighborhood of one or the other of the two ori- 
entations, perhaps with occasional reversals of the flow 
direction. 

A long-term goal of the approach that we have taken 
would be to apply these results to natural systems, such 
as flow in Earth's outer core or the solar convection zone 
which are spherically symmetric shells with a Coriolis 
force. These flow structures are much more complicated 
than a single LSC roll, but while the asymmetries align 
these circulations, we have shown that when these asym- 
metries are small they do not affect the LSC strength and 
should not affect the frequency of cessations. However, 
further experimentation should be carried out to estab- 
lish whether stronger asymmetries that align the LSC 
can also affect its stability to cessations. 



VIII. APPENDIX: COMPARISON OF 
DIFFUSIVITIES 

The basic model Eqs.[T] and [2] involves the two stochas- 
tic forces fs {t) and fg (t) that lead to diffusivities Dg and 
Dg of S and 9. Here we compare these diffusivities. One 
might hypothesize that the thermal fluctuations in the 
system are isotropic, so the fluctuations would have the 
same strength both in the direction of the LSC velocity 
and in the azimuthal motion. From Eq. [5]the variance of 
temperature fluctuations is given by 

erf, ^ ((^5)2 ^Qg2(g _ _^ {dOo)^S^ sin2(6l - 6'o)X47) 
« [(dS)^ + id9ofP]/2 . (48) 

If the S and Oq terms contribute equally to the total fluc- 
tuations, then equipartition gives (dS)'^ = S^{d9o)^. Us- 
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simplifying to the /3 = case where 5 = Sq leads to 



FIG. 12: The dimensionless ratio between diffusivities 
Ds/{DffS^T?) for (3 — 0. Solid symbols; data from the 
L = 24.76 cm sample. Open symbols: data from a similar 
sample with L = 50.61 cm, referred to as the large sample in 
Ref. [2^. The ratio is expected to be 1 if fluctuations have 
the same strength in both the (j) and 9 coordinates. 



ing the difFusivity relationships found in Ref. [23| and 



(49) 



We plot the 1. h. s. of this equation in Fig. [12] using the 
experimental values from Ref. psj for /3 = 0. The ratio 
is slightly less than 1 with little variation over our ex- 
perimental range of i?, so by this measure the thermal 
fluctuations in the direction of the LSC are smaller than 
in the azimuthal direction, if only slightly. For the pur- 
poses of the model, the fluctuations can be considered 
nearly isotropic, so use of Eq. [49] reduces the number of 
free parameters of the model by 1. 



IX. ACKNOWLEDGMENTS 

This work was supported by Grant DMR07-02111 of 
the US National Science Foundation. We thank Chris 
Henley for his insightful analysis of our data. 



[1] E. D. Siggia, High Rayleigh number convection, Aiirm. 

Rev. Fluid Mech. 26, 137 (1994). 
[2] L. P. Kadanoff, Turbulent heat flow: Structures and scal- 
ing, Phys. Today 54, 34 (2001). 
[3] G. Ahlers, S. Grossmann, and D. Lohse, Hochprazision 

im Kochtopf: Neues zur turbulenten Konvektion, Physik 

Journal 1, 31 (2002). 
[4] F. Heslot, B. Castaing, and A. Libchaber, Transition to 

turbulence in helium gas, Phys. Rev. A 36, 5870 (1987). 
[5] M. Sano, X. Z. Wu, and A. Libchaber, Turbulence in 

helium-gas free convection, Phys. Rev. A 40, 6421 (1989). 
[6] B. Castaing, G. Gunaratne, F. Heslot, L. Kadanoff, A. 

Libchaber, S. Thomae, X. Z. Wu, S. Zaleski, and G. 

Zanetti, Scaling of hard thermal turbulence in Rayleigh- 

Benard convection, J. Fluid Mech. 204, 1 (1989). 
[7] S. Ciliberto, S. Cioni, and C. Laroche, Large-scale flow 

properties of turbulent thermal convection, Phys. Rev. E 

54, R5901 (1996). 
[8] T. Takeshita, T. Segawa, J. A. Glazier, and M. Sano, 

Thermal turbulence in mercury, Phys. Rev. Lett. 76, 

1465 (1996). 

[9] S. Cioni, S. Ciliberto, and J. Sommeria, Strongly turbu- 
lent Rayleigh- Benard convection in mercury: comparison 
with results at moderate Prandtl number, J. Fluid Mech. 
335, 111 (1997). 

[10] X. L. Qiu and P. Tong, Large-scale coherent rotation and 
oscillation in turbulent thermal convection, Phys. Rev. E 
61, R6075 (2000). 

[11] X. L. Qiu and P. Tong, Large scale velocity structures in 
turbulent thermal convection, Phys. Rev. E 64, 036304 
(2001). 

[12] X. L. Qiu and P. Tong, Onset of coherent oscillations in 
turbulent Rayleigh- Benard convection, Phys. Rev. Lett 
87, 094501 (2001). 

[13] J. J. Niemela, L. Skrbek, K. R. Sreenivasan, and R. J. 



Donnelly, The wind in confined thermal turbulence, J. 
Fluid Mech. 449, 169 (2001). 
[14] X. L. Qiu and P. Tong, Temperature oscillations in tur- 
bulent Rayleigh- Benard convection, Phys. Rev. E 66, 
026308 (2002). 

[15] X. L. Qiu, X. D. Shang, P. Tong, and K.-Q. Xia, Veloc- 
ity oscillations in turbulent Rayleigh- Benard convection, 
Phys. Fluids. 16, 412 (2004). 

[16] D. Funfschilling and G. Ahlers, Plume motion and large 
scale circulation in a cylindrical Rayleigh- Benard cell, 
Phys. Rev. Lett. 92, 194502 (2004). 

[17] C. Sun, K. Q. Xia, and P. Tong, Three-dimensional flow 
structures and dynamics of turbulent thermal convection 
m a cylindrical cell, Phys. Rev. E 72, 026302 (2005). 

[18] Y. Tsuji, T. Mizuno, T. Mashiko, and M. Sano, Mean 
Wind in Convective Turbulence of Mercury, Phys. Rev. 
Lett. 94, 034501 (2005). 

[19] C. Sun, H. D. Xi, and K. Q. Xia, Azimuthal symmetry, 
flow dynamics, and heat transport in turbulent thermal 
convection in a cylinder with an aspect ratio of 0.5, Phys. 
Rev. Lett. 95, 074502 (2005). 

[20] H. D. Xi, Q. Zhou, and K. Q. Xia, Azimuthal motion 
of the mean wind in turbulent thermal convestion, Phys. 
Rev. E 73, 056312 (2006). 

[21] E. Brown and G. Ahlers, Rotations and cessations of the 
large-scale circulation in turbulent Rayleigh- Benard con- 
vection, J. Fluid Mech. 568, 351 (2006). 

[22] E. Brown and G. Ahlers, Effect of the Earth's Coriolis 
force on turbulent Rayleigh-Benard convection in the lab- 
oratory, Phys. Fluids 18, 125108 (2006). 

[23] E. Brown, A. Nikolaenko, and G. Ahlers, Reorientation of 
the large-scale circulation in turbulent Rayleigh- Benard 
convection, Phys. Rev. Lett 95, 084503 (2005). 

[24] S. Q. Zhou, C. Sun, and K. Q. Xia, Measured oscilla- 
tions of the velocity and temperature fields in turbulent 



16 



Rayleigh-Benard convection in a rectangular cell, Phys. 

Rev. E 76, 036301 (2007). [39] 

C. Resagk, R. du Puits, A. Thess, F. V. Dolzhansky, S. 
Grossmann, F. Fontenele Araujo, and D. Lohse, Oscilla- 
tions of the large scale wind in turbulent thermal convec- 
tion, Phys. Fluids 18, 095105 (2006). [40] 

D. Funfschilling, E. Brown, and G. Ahlers, Azimuthal 
oscillations of the large-scale circulation in turbulent 
Rayleigh-Benard convection, J. Fluid Mech 607, 119 
(2008). [41] 

E. Brown and G. Ahlers, Large-scale circulation model of 
turbulent Rayleigh-Benard convection, Phys. Rev. Lett. 

98, 134501 (2007). [42] 
E. Brown and G. Ahlers, A model of diffusion in a poten- 
tial well for the dynamics of the large-scale circulation in 
turbulent Rayleigh-Benard convection, Phys. Fluids 20, [43] 
075101 (2008). 

E. Brown, D. Funfschilling, and G. Ahlers, Anomalous 
Reynolds-number scaling in turbulent Rayleigh-Benard 
convection, J. Stat. Mech. P10005 (2007). [44] 
G. Glatzmaier, L. H. R. Coe, and P. Roberts, The role 
of the Earth's mantle in controlling the frequency of geo- 
magnetic reversals, Nature(London) 401, 885 (1999). [45] 
R. Howard and B. LaBonte, The Sun is observed to be a 
torsional oscillator with a period of 11 years, Astrophys. 
J. 239, L33 (1980). 

E. van Doorn, B. Dhruva, K. R. Sreenivasan, and V. [46] 
Cassella, Statistics of wind direction and its increments, 

Phys. Fluids 12, 1529 (2000). 

G. Ahlers, E. Brown, and A. Nikolaenko, The search for 
slow transients, and the effect of imperfect vertical align- 
ment, in turbulent Rayleigh-Benard convection, J. Fluid [47] 
Mech. 557, 347 (2006). 

F. Chilla, M. Rastello, S. Chaumat, and B. Castaing, 

Long relaxation times and tilt sensitivity in Rayleigh- [48] 

Benard turbulence. Euro. Phys. J. B 40, 223 (2004). 

J. J. Niemela and K. R. Sreenivasan, Confined turbulent 

convection, J. Fluid Mech. 481, 355 (2003). [49] 

F. Fontenele Araujo, S. Grossmann, and D. Lohse, Wind 

reversals in turbulent Rayleigh-Benard convection, Phys. 

Rev. Lett. 95, 084502 (2005). [50] 

R. Benzi, Flow reversal in a simple dynamical model of 

turbulence, Phys. Rev. Lett. 95, 024502 (2005). 

M. van Reeuwijk, H. Jonkers, and K. Hanjalic, Wind and 

boundary layers in Rayleigh- Bnard convection. L Analy- 



sis and modeling, Phys. Rev. E 77, 036311 (2008). 
M. van Reeuwijk, H. Jonkers, and K. Hanjalic, Wind 
and boundary layers in Rayleigh-Bnard convection. IL 
Boundary layer character and scaling, Phys. Rev. E 77, 
036312 (2008). 

E. Brown, A. Nikolaenko, D. Funfschilling, and G. 
Ahlers, Heat transport in turbulent Rayleigh-Benard con- 
vection: Effect of finite top- and bottom-plate conductiv- 
ities, Phys. Fluids 17, 075108 (2005). 
H. D. Xi and K. Q. Xia, Cessations and reversals of 
the large-scale circulation in turbulent thermal convec- 
tion, Phys. Rev. E 75, 066307 (2007). 
H. D. Xi and K. Q. Xia, Flow mode transitions in tur- 
bulent thermal convection, submitted to Phys. Fluids 
(2007). 

E. Brown and G. Ahlers, Temperature gradients, and 
search for non-Boussinesq effects, in the interior of tur- 
bulent Rayleigh-Benard convection, Europhys. Lett. 80, 
14001 (2007). 

S. Grossmann and D. Lohse, Prandtl and Rayleigh num- 
ber dependence of the Reynolds number in turbulent ther- 
mal convection, Phys. Rev. E 66, 016305 (2002). 
A. Nikolaenko, E. Brown, D. Funfschilling, and G. 
Ahlers, Heat transport by turbulent Rayleigh-Benard con- 
vection in cylindrical cells with aspect ratio one and less, 
J. Fluid Mech. 523, 251 (2005). 

In Ref. [2^ the potential Eq. [13] was found to be incon- 
sistent with the data for p{5) at small 5 < 0.55o- Thus an 
empirical potential was introduced for which, consistent 
with the data, dV/dS = —BDsSo for small 5 with B = 12 
for = 1.1 X 10"'. 

H. A. Kramers, Brownian motion in a field of force and 
the diffusion model of chemical reactions, Physica (Ams- 
terdam) 7, 284 (1940). 

M. Gitterman, The Noisy Oscillator, The First Hundred 
Years, From Einstein Until Now (World Scientific, Sin- 
gapore, ADDRESS, 2005). 

Z. A. Daya and R. E. Ecke, Does turbulent convection feel 
the shape of the container?, Phys. Rev. Lett. 87, 184501 
(2001). 

R. Krishnamurti and L. N. Howard, Large scale flow gen- 
eration in turbulent convection, Proc. Natl. Acad. Sci. 78, 
1981 (1981). 



